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LIGHT  SCATTERING  FROM  ROUGH  SURFACES  ! 
FINAL  REPORT 
EXECUTIVE  SUMMARY 


The  particular  aspect  of  light  scattering  investigated  in  this  project  was  the  angular  correlation 
of  speckle  patterns  produced  by  scattering  from  rough  surfaces.  The  Report  consists  primarily 
of  a  draft  copy  of  the  Thesis  of  R  J  Syratt,  the  student  employed  on  the  Contract 

A  randomly  rough  surface  illuminated  by  a  parallel  laser  beam  at  some  angle  0  produces  a 
speckle  pattern  in  space.  If  the  angle  of  illumination  is  changed  to  0  +  50  then  the  speckle 
pattern  changes  in  two  ways:  (i)  it  translates  in  angle  and  (ii)  it  changes  its  \  m  in  a  manner 
that  depends  on  the  detailed  surface  characteristics  and  the  angles  involved.  We  have  made 
experimental  measurements  of  the  angular  correlation  of  speckle  for  two-dimensional  Gaussian 
surfaces  of  low  and  high  root-mcan-squarc  slope.  In  comparing  with  theory  we  have 
distinguished  between  surfaces  for  which  single  scattering  is  valid  (very  low  slope)  and  those 
for  which  Kirchhoff  theory  is  valid  (larger  rms  slopes).  A  numerical  ray-tracing  approach  has 
been  implemented  for  two-dimensional  surfaces  that  includes  shadowing  and  polarisation  in  a 
limited  way.  A  conceptual,  heuristic  model  is  developed  that  explains  the  essential  physical 
features. 

In  addition  to  the  main  subject  of  this  study,  we  also  investigated  a  number  of  subsidiary  topics 
as  detailed  in  the  attached  publications. 
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Multiple  scattering  from  rough  dielectric  and  metal 
surfaces  using  the  Kirchhoff  approximation 

N,  C.  BRUCE  and  J.  C.  DAINTY 
The  Blackett  Laboratory,  Imperial  College, 

London  SW7  2BZ,  England 

(Received  11  Jataury  1991) 

Abstract.  The  Kirchhoff  double-scatter  method  for  ralculating  the  intensity 
distribution  scattered  from  a  rough  surface  is  extended  to  dielectric  and  metal 
surface  materials.  The  material  properties  are  contained  in  the  Fresnel  reflection 
coefficients  only.  It  is  shown  that  the  results  agree  well  with  calculations  using  the 
exact  method  for  a  surface  of  Gaussian  statistics  with  standard  deviation  of  height 
u=  1-931  and  1/e  correlation  length  x = 5-02L 


L  Introduction 

In  a  previous  paper  [1]  it  was  shown  that  using  the  Kirchhoff  approximation 
(KA)  and  including  shadowring,  it  was  possible  to  derive  equations  for  the  single-  and 
double-scatter  contributions  to  the  intensity  distribution  scattered  from  a  randomly 
rough  surface.  In  [1],  the  simplest  situation  was  considered — that  of  a  perfectly 
conducting  surface.  As  expected,  the  double-scatter  term  showed  the  enhanced 
backscatter  effect  [2-fl,  with,  a  peak  approximately  equal  to  twice  the  background  in 
the  backscatter  direction.  This  agrees  with  the  simple  ray  picture  of  the  scattering 
process  for  which  rays  and  their  time-reversed  partners  add  coherently  in  the 
backscatter  direction  and  incoherently  in  other  directions. 

A  comparison  of  the  perfect  conductor  calculations  with  experimental  distri¬ 
butions  of  the  scattered  light  from  a  go'd-coated  surface  with  Gaussian  statistics 
showed  good  agreement  at  low  angles  of  incidence  (up  to  about  20°  from  normal)  but 
poor  agreement  at  higher  angles  [1].  The  same  trend  is  seen  in  comparisons  of  exact 
calculations  using  die  extinction  theorem  with  experimental  values  [4],  although  for 
perfectly  conducting  surfaces  the  calculated  intensities  for  die  extinction  and 
Kirchhoff  double-scatter  methods  agree  very  well  [5].  This  suggests  that  the 
Kiichhpff  method  has  included  the  main  physical  processes  that  produce  enhanced 
backscatter. 

In  this  paper  the  results  of  the  method  for  the  more  general  cases  of  scattering 
from  a  dielectric  and  a  general  conductor  (metal)  are  reported.  In  section  2  the  theory 
is  briefly  outlined  and  in  section  3  the  results  are  shown  and  compared  with  exact 
calculations  and  experimental  results. 

2.  Theory 

The  starting  point  is  the  two-dimensional  Helmholtz  integral  equation  for  a 
scattered  field  (the  surface  profile  is  constant  along  the  y-direction): 
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where  Et(x,  s')  is  the  total  field  at  the  point  (*',  *0.  r = [(* —  *')2 + (* — *')2] 1/2 > 

is  the  zeroth  order  Hankel  function  of  the  first  kind  and  dr'  is  an  element  of  the 

surface. 

The  KA  approximates  the  total  field  on  the  surface  as  the  sum  of  the  incident  field 

plus  the  reflected,  field  at  that  point: 

E,(x.z)=(l  +R)Ei(x,z),  (2) 

dEfaz) 


dn 


*.i(l-/i)k,  •«£,(*.*)» 


(3) 


where  R  is  the  planar  reflection  coefficient  at  that  point,  depending  on  the  local 
incidence  angle;  Et(xy  g)  is  the  field  incident  at  that  point;  kj  is  the  incident  wave- 
vector;  and  n  is  the  outward  normal  to  the  surface  at  the  point  ( x , Since  R  is  the 
planar  reflection  coefficient,  we  immediately  have  a  condition  for  the  validity  of  the 
approximation — the  surface  must  be  locally  flat  or,  as  is  usually  written,  the  radius  of 
curvature  of  the  surface  must  be  large  compared  with  the  wavelength. 

Substituting  (2)  and  (3)  into  (1)  and  performing  a  small  amount  of  mathematics, 
following  Beckmann  and  Spizzichino  [6],  the  standard  single-scatter  solution  is 
obtained; 


where  the  reflectivity  is  a  funedor  of  x,  since  die  local  angle  of  incidence  varies  along 
die  surface".  In  equation  (4),  fl,  :s  the  incidence  angle,  8  is  the  angle  of  scatter 
(measured  positive  in  the  opposite  direction  to  the  incident  angle)  and  <p  is  a  phase 
factor  depending  only  on  (x,  sr)  and  not  on  (x\  s').  To  ohtain  this  expression  we  have 
assumed  an  incident  plane  wave  of  unit  amplitude.  Shadowing  effects  must 
somehow  be  included  in  this  method .»  have  a  physically  realistic  result.  This  is  done 
by  multiplying  die  integrand  in  (4)  by  incidence  and  scatter  shadow’  functions 

*  f 1  if  {*',*')  is  illuminated, 

(0  i-- (x', af)  is  not  illuminated, 

fl  (*'>*')  is  visible, 

*  '  *  [0  if  (*',*')  is  not  visible. 

These  represent  geomedical  or  straight  line  shadow  functions  and  so  are  an 
approximation  to  die  true  effect  of  shadowing. 

The  derivation  of  the  double-scatter  contribution  involves  two  terms.  The  first  is 
the  field  scattered  from  one  point  on  the  surface  to  another  point,  also  on  the  surface, 

£.<>„.,>=  -i(<l  +*.)(«.  HVW.,) 

-(1  -RJiMmi  sin#i+cos0|)/iVl(fcr,,)j£|(*i,*i),  (S) 

where  subscript  1  represents  the  first  point  and  subscript  2  the  second.  The  normal 
derivative  of  the  Hankel  function  is  given  by 
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The  second  term  is  the  field  scattered  into  the  space  above  the  surface  from  every 
combination  of  first  and  second  points: 

£?(#)= ~ —■  Jf  J  S\ 2^0  +  sin 0- cos 0)E,(x2, z2) 


■fi(l -R2)  ^^'^-^expih^Radxi  dxj, 


(7) 


where  we  have  included  another  shadowing  function  todescribe  shadowing  between 
points  on  the  surface: 


^  _  fl  if  (x2,z2)  is  visible  1 
12  (0  if  (x2,Z2)  is  not  visi 


from  (*„*,), 
visible  from  (xl,zl), 


and  substituted 
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In  equation  (7)  we  have  used  equation  (2)  a  second  time  to  write  that  the  total  field  at 
the  second  point  is  the  field  from  the  first  point  plus  its  reflection  (J?2).  This  can  be 
seen  to  be  reasonable  if  it  is  noted  that  the  condition  on  the  radius  of  curvature 
depends  only  on  the  surface  profile  and  so  will  be  true  for  all  subsequent  interactions 
if  it  is  true  for  the  first.  The  validity  of  the  parameters  used  here  was  discussed  in  the 
previous  paper.  This  method  is  the  same  as  the  iterative  Kirchhoff  solution  {7, 8]  but 
with  the  inclusion  of  the  shadow  functions  giving  the  effect  of  shadowing  explicitly 
rather  than  implicitly  as  in  the  iterative  method  [8], 

Equations  (5)  and  (7)  together  form  the  double-scatter  contribution.  The 
simplest  case  to  consider  is  when  die  material  is  perfectly  conducting.  Then  the 
reflection  coefficients  R  are  either  1  or  —1  for  p  (TM)  and  s  (TE)  polarizations, 
respectively,  and  we  are  left  with  only  one  term  as  already  discussed.  For  the  more 
general  case  we  end  up  with  the  sum  of  four  terms.  In  all  the  cases  considered,  the 
normal  derivative  dE,(x2,  x 2)ldn2  is  approximated  in  the  following  way: 

dEn(x2, zz)  ^  1  /(1  ,  ^  J  fe(z2— 

dn2  4iV  i;V  1  ri2  )  d*>2 

-(1  -jR1)iMw1  sinfli+cosfl,)  E-,(xu  zt).  (8) 

In  (8)  we  need  to  use 


aff?>(*T12)  f  A(X2- Xt)  ,  _  rrWIW  ,  v, 

p..  —J  nix  +rt2i  hrltfo  (*fi2)— Hi  (^la)]* 

2  l.  ^12  t^X2  J2 


The  material  dependence  of  these  equations  is  contained  solely  in  the  Fresnel 
reflection  coefficients  R.  For  scattering  from  a  dielectric  surface  of  refractive  index  n 
in  air  («=  1),  the  reflection  coefficients  for  p  and  s  polarizations  are  [9] 


noosflj— 

*r= - 

n  cos  6 1 + 
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cos  et+nj(l-^j 


To  find  the  scattered  field  we  substitute  either  (9)  or  (10),  depending  on  the 
polarization,  into  (S)  and  (7).  This  case  of  the  perfect  dielectric  is  reasonably 
straightforward.  A  more  complicated  situation  occurs  when  the  case  of  light 
scattered  from  an  interface  between  a  dielectric  and  air  is  considered  when  the  light  is 
incident  from  the  dielectric  side  of  the  boundary.  In  this  case  the  Fresnel  coefficients 
are  [9] 

p  c<>s  Oi-njl-n1  sin2  6, 

^  cos  Oj  +  n^/l  —it2  sin2#;’ 


Then,  when 


n  cos  —  y/l  —  n2  sin2  6t 
n  cos  0, + ^/l  —  n2  sin2  0; 

n2  sin2  0j  >  1 , 


the  reflectivities  are  complex.  This  occurs  when  the  critical  incidence  angle  0C  is 
reached  and  we  lave  total  internal  reflection 


The  final  case  considered  is  that  of  die  general  conductor  or  metal  which  involves 
a  complex  refractive  index  «-*«4-iK.  In  this  case  the  reflectivities  are  [9] 

(n2  cos2  0,-q2  cos2  y)  4-  (k2  cos2  fft —q2  sin2  y) 

^  (ncosOt+qcosy'f+iKcosOi+q&iny)2 

t  2.  ((cgcosgicosy— ngsinycosg;)  (13) 

(a  cos  0; 4- q  cos  y)2 + (K  cos  0j + q  sin  y)2  * 

_  cos2  0; — (n^cos  y—Kq  sin  y)2 — (icq  cos  y+ttq  sin  y)2 
(cos  0j  4-  (nq  cos  y — xq  sin  y))2 + {Kq  cos  y  +  nq  sin  y)2 

_2_ _ cos  9j(Kq  cosy+nq  sin  y) _ 

(tosdl+(nqcosy—Kq$my))2+(Kqcosy+nqsmy)2’ 

where 

«  (,  (h2  —  *2)  .  (  2mc  .  -„\2  .... 

»  Tm"  *■)  8‘)  ■  ,,5) 

(n2-x2)  . 

tan  2y =- — nf>\ 


(«2  +  k2)2 


It  should  be  noted  that  the  contribution  from  light  paths  that  traverse  inside  parts 
of  the  surface  are  not  taken  into  account.  It  is  believed  that  such  terms  are  small  due 
to  the  requirement  for  large  angles  of  deviation  of  refraction  to  diiect  tight  back  into 
the  space  above  the  surface. 
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3.  Results 

The  Kirchhoff  method  equations  have  been  discretised  in  the  usual  way 
(assuming  that  values  of  parameters  remain  constant  over  a  suitably  small  range  of 
variables)  and  programmed  into  a  Sun  4  computer.  A  surface  of  length  302  was 
discretized  into  200  segments  giving  a  segment  length  of  01 52.  Computation  time? 
were  7  minutes  per  frame  for  die  dielectric  cases  and  10  minutes  per  frame  for  the 
metal  with  approximately  600  frames  averaged  for  each  graph.  Computational 
results  are  shown  in  figure  1  with  experimental  results  [10]  shown  in  figure  2  for 
comparison.  The  surface  parameters  are:  standard  deviation  of  height  o  =  1-18  pm, 
correlation  length  t  =  2-97  pm,  with  a  wavelength  of  0  633  pm  (He-Ne  red)  and  a 
refractive  index  of  n  —  1  -41 .  The  scattered  energy  as  a  function  of  incidence  angle  for 
both  methods  is  given  in  table  1. 

The  Kirchhoff  and  exact  methods  agree  quite  well  both  for  scattered  energy  and 
the  intensity  distribution  up  to  —30°  incidence.  In  both  the  s  and  p  cases  the 
Kirchhoff  scattered  intensity  distribution  is  rising  at  scatter  angles  near  +  90°  and  for 
p-polarization  the  energy  scattered  is  twice  its  value  for  the  extinction  case.  The 
cause  of  this  may  be  that  the  effect  of  light  paths  that  pass  through  the  material  are 
becoming  important,  or  the  approximation  used  for  the  shadow  functions,  i.e. 
straight  line  or  geometric  shadow  functions,  is  less  valid  when  incidence  shadowing 
starts  to  have  an  effect.  The  second  reason  is  supported  by  the  fact  that  for  the 
perfectly  conducting  case  the  unitarity  was  from  6  to  9%  greater  than  unity.  There  is 
reasonable  agreement  with  the  shape  of  the  experimental  curves,  although  note  that 
the  experimental  values  have  not  been  normalized.  However,  from  the  Kirchhoff 
results  it  can  be  seen  that  the  double-scattered  energy  is  an  order  of  magnitude  down 
on  the  single  scatter  for  s-polarization  (TE)  and  several  orders  of  magnitude  down 
for  p-polarization  (TM).  This  means  that  scattering  from  this  dielectric  surface  is 
mainly  a  single-scatter  effect,  so  that  even  though  enhanced  backscatter  can  be  seen 
in  the  double-scatter  curves,  the  effect  on  the  total  scatter  curve  is  very  small..  This  is 
probably  because  only  a  small  part  of  the  energy  is  reflected  at  each  interaction  point 
on  the  surface. 

Assuming  that  in  this  single-scatter  regime  the  dominant  term  from  each  part 
of  the  surface  is  the  specular  term  it  can  be  shown  that  the  minimum  in  the 
p-polarization  case  is  directly  related  to  the  Brewster  angle.  Figure  3  shows  the 
situation  where  one  of  the  discretized  surface  facets  is  illuminated  at  the  Brewster 
angle  and  so  will  not  scatter  p-polarized  light.  Then  6„  the  scatter  angle,  is  given  by 

—  0, — Q\  4-  2  0B, 

since  8t  is  negative.  For  n— 1-41,  the  Brewster  angle  is 

0B  =  tan~‘(l-41)ss55°. 

I  bus  at  0i  —  —  30°,  0,=  —80°  and  at  0(  =  —60°,  0,=  —50°,  which  agrees  reasonably 
well  with  the  observed  ingle  at  which  the  minimum  occurs.  Another  feature  of  the 
curves  which  can  be  explained  using  the  single-scatter  model  is  that  the  s  and  p 
curves  have  the  same  value  at  backscatter  for  all  incidence  angles.  If  the  specular 
term  is  dominant  then  the  backscatter  results  from  parts  of  the  surface  which  are 
normal  to  the  incident  direction.  The  Fresnel  coefficients  for  normal  reflections  have 
the  same  modulus  but  different  sign  for  s  and  p,  so  the  scattered  intensities  are  the 
Same  in  the  backscatter  direction. 
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Figu>  c  2.  Experimental  curves  for  scattering  from  a  dielectric  as  in  figure  1 :  s  Cupper  curve) 
and  p  (lower  cut  ve)  polarizations. 


Table  i.  Scattered  energy  as  a  function  of  incidence  angle  (dielectric). 
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Figure  3.  Geometry  tor  scattering  at  Brewster  angle. 


The  results  for  the  inverse  dielectric  case  for  the  Kirchhoff  method  with 
s-polarization  and  the  same  material  as  above  («=  1-41)  are  shown  in  figuie  4.  The 
scattered  energy  is  shown  in  table  2. 

In  this  case,  the  double  scatter  forms  a  much  larger  part  of  the  total  scattered 
energy  so  that  a  strong  enhanced  backscatter  peak  is  expected.  Indeed,  this  is  seen  at 
low  angles  of  incidence  but  the  peak  dies  away  quickly  so  that  by  —  20°  it  has  almost 
disappeared  completely.  The  steep  slope  of  the  double-scatter  curve  between  0°  and 
20°  gives  this  effect. 

The  peak  in  the  single-scatter  curves  at  large  negative  angles  can  be  explained  in 
the  same  way  as  the  Brewster  angle  effect  above.  The  scatter  angle  corresponding  to  a 
local  angle  of  incidence  equal  to  the  critical  angle  is  given  by 

-$,-9,  +  20„ 

where  the  critical  angle  6C  is  given  by 


giving  6%  =  —  70°  for  8-,  =  —  20°..  This  angle  should  give  the  start  of  the  region  of  total 
internal  reflection.  Shadowing  counteracts  this  effect  by  blocking  light  scattered  at 
high  angles,  causing  the  reflectivities  to  come  down  again  toward  —  90°. 

The  final  case  considered  is  that  of  scattering  from  a  metal,  in  particular  gold,  at 
2=0-633  pm,  with  the  refractive  index  n  =  0167-Fi3149  [11],  Previous  comparisons 
of  experiment  and  theory  have  used  calculations  from  perfect  conductors,  so  it  would 
be  useful  to  find  the  effect  of  finite  conductivity  on  the  scattered  intensity.  The 
scattered  energy  for  two  angles  of  incidence  is  given  in  table  3. 
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Table  3.  Scattered  energy  for  two  different  angles  of  incidence  (gold). 


Incidence 

angle 

Kirchhoff 

single 

Kirchhoff 

double 

Kirchhoff 

total 

s-polarisation  ( TE ) 

0° 

0-811 

0-105 

0-902 

-40° 

0-847 

0-075 

0  872 

p-polarization  (TM) 

0° 

0-783 

0-112 

0-882 

—  40° 

0-829 

0-069 

0-849 

p  pol.  s  po!. 


Figure  5.  As  figure  1.  for  a  metal  surface  at  0°  (top)  and  —40°  incidence.  Kirchhoff  total 
(solid  line),  single  (oooo)  and  double  (  +  +  +  +)  renormalized  to  have  unit  area  under 
the  total  curve.  The  dotted  line  is  the  Kirchhoff  total  for  the  same  surface  as  a  perfect 
conductor. 
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These  values  of  scattered  energy  could  be  used  as  a  check  on  the  parameters  of  the 
surface  material,  allowing  a  check  on  whether  there  is  an  unexpected  effect  present  in 
the  scattering  process,  for  example,  if  the  metal  coating  is  too  thin  or  if  there  is  some 
granulation  of  the  surface  material.  Here,  to  compare  the  perfect  conductor  and 
metal  results,  the  two  cases  are  normalized  to  the  same  area  under  the  graph  (as  has 
been  done  on  previous  experimental  results)  and  compared  directly.  Figure  5  shows 
the  results.  There  is  very  little  difference  in  the  metal  and  perfectly  conducting 
curves,  showing  that  finite  conductivity  has  very  little  effect  on  the  scattered 
intensity  distribution. 


4.  Conclusions 

Using  the  Kirchhoff  double-scatter  contribution  it  has  been  shown  that 
scattering  from  a  particular  dielectic  surface  is  mainly  a  single-scatter  effect  when  the 
light  is  incident  from  the  air  side,  but  when  the  incidence  is  from  the  dielectric  side 
much  more  light  is  double  scattered.  We  believe  that  the  reflectivities  at  each 
interaction  are  very  low  in  the  first  case  and  total  internal  reflection  gives  relatively 
more  double  scattering  in  the  second..  The  results  agree  well  with  exact  calculations 
using  the  extinction  theorem,  except  at  higher  angles  of  incidence  where  it  is  thought 
that  the  approximation  used  for  shadowing  becomes  less  valid.  Scattering  from  a 
metal  surface  was  compared  with  scattering  from  a  perfect  conductor  and  found  to 
agree  well.  This  shows  that  the  effect  of  finite  conductivity,  at  least  for  very  good 
conductors,  is  purely  to  reduce  the  scattered  energy  and  not  to  alter  the  scattered 
intensity  distribution. 
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Abstract  Measurements  are  presented  of  the  angular  distribution  of  four  wavelength'  of 
light  scattered  by  a  one-dimensional  random  rough  surface,  whose  probability  density 
function  is  Gaussian  with  a  standard  deviation  a  =  1.22  +  0.02  tun  and  whose  lateral 
correlation  function  is  also  Gaussian  with  1/e  width  t  =  3.17  +  007  pm  The  wavelengths 
used  are  0.63, 1.15, 3.39  and  10.6  pm.  The  surface  is  used  in  two  forms:  coated  with  gold 
and  as  an  almost  lossless  dielectric.  The  results  are  compared  to  those  predicted  by  a 
double  scattering  form  of  the  Kirchhoff  formulation.  Agreement  is  good  at  small  angles  of 
incidence  but  less  good  at  larger  angles  of  incidence. 


1.  Introduction 

The  experimental  observation  of  enhanced  backscaltering  from  random  rough  surfaces 
of  large  root-mean-square  slope,  first  reported  by  Mendez  and  O’Donnell  [  1,2],  has 
stimulated  a  re-examination  of  the  problem  of  light  scattering  in  the  past  few  years.  The 
main  progress  to  date  has  been  the  development  of  ‘exact’  numerical  codes  for  the 
solution  of  Maxwell's  equations  from  a  one-dimensional  surface  illuminated  with  cither  s 
(i.e.  te)  polarization  or  p  (ie.  tm)  polarization  [3-7],  With  s  polarization,  the  electric 
vector  is  parallel  to  the  grooves,  whereas  with  p  polarization  it  is  perpendicular  to  the 
grooves,  as  in  figure  1  (this  figure  also  shows  the  sign  convention  used  for  the  incident 
and  scattering  angles).  An  important  feature  of  the  work  of  Meudez  and  O’Donnell 
was  that  the  surfaces  were  relatively  well  characterized,  with  Gaussian  statistics  for  the 
surface  height  and  a  single-scale  Gaussian  correlation  function.  Since  the  statistics  of 
the  surface  were  known,  a  critical  comparison  between  experiment  and  theory  could  be 
made  with  confidence. 

The  shape  of  the  scattering  cross  section  curves  with  angle  of  observation  for 
high-sloped  surfaces  is  quite  different  from  that  for  simple  low-sloped  ones  and  early 
numerical  results  [3]  were  encouraging  since  they  were  in  fairly  good  agreement  with 
the  experimental  ones  particularly  at  near-normal  incidence.  In  order  to  carry  out  a 
more  critical  comparison  between  real  experiments  and  numerical  ones,  it  is  important 
that  the  surface  is  well  characterized  and  also  helpful  if  a  range  of  wavelengths  are  used. 
The  results  presented  here  are  intended  to  supplement  those  already  reported  [8-10] 
with  the  aim  of  providing  a  reliable  body  of  experimental  data  for  comparison  with 
numerical  work  and  analytical  theory.  The  surface  used  is  one-dimensional,  for  two 
reasons:  firstly,  it  can  be  characterized  much  more  accurately  than  a  two-dimensional 
one,  since  a  sharp  chisel-shaped  stylus  can  be  used  in  a  surface  profilometer;  secondly, 
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Figure  1.  Polarization  and  angle  notation  tor  in-planc  scattering  from  a  onc-dimensional 
rough  surface 


‘exact’  numerical  calculations  of  the  light  scattering  are  only  feasible  at  the  present  time 
for  the  one-dimensional  case. 

When  comparing  experimental  measurements  of  light  scattering  with  numerical 
computations,  it  is  heipful  if  the  numerical  results  provide  some  physical  insight  to  the 
scattering  process.  For  example,  it  is  believed  from  the  experiments  that  the  mechanism 
giving  rise  tc  the  enhanced  backscatter  peak  is  multiple  scatrering;  however,  numerical 
calculations  based  on  exact  theory  do  not  separate  the  single  and  multiple  scatter  terms, 
and  therefore  do  not  provide  the  insight  that  is  desirable  (however,  iterative  solutions 
do  separate  the  single  and  multiple  scatter  terms).  For  this  reason,  we  have  written 
numerical  code  based  on  a  multiple  (double  and  triple)  scattering  extension  of  the 
Kirchhoff  boundary  condition,  including  the  effects  of  shadowing  (see  [11]  for  details 
and  further  references).  In  section  3  of  this  paper  we  compare  the  results  of  this  code 
with  the  experimental  results  and  ‘exact’  numerical  code. 


2.  Experimental  results 

Master  surfaces  are  produced  by  exposing  a  thick  layer  of  photoresist  ( as  12  /an  of 
Shipley  S1400-37)  to  several  statistically  independent  laser  speckle  patterns.  Two  versions 
of  the  surface  were  prepared  using  a  replication  technique  [8],  one  being  coated  with 
a 1000  A  of  gold  and  the  other  being  an  almost  lossless  dielectric  of  refractive  index 
n  —  1.41  (at  k  —  0.63  /tm).  Figure  2  shows  the  probability  histogram  of  surface  height 
and  surface  autocorrelation  function,  as  measured  by  a  Talystep  profilometer  whose 
stylus  is  a  pyramid  of  70°  apex  angle  truncated  by  a  flat  region  of  a  0.5  /an.  Both  are 
good  fits  to  Gaussian  functions,  with  the  root-mean-square  height  a  =  1.22  ±  0.02  /an 
and  1/e  correlation  lengtn  t  =  3.17  ±  0.07  /an.  The  angular  distribution  of  the  scattered 
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Figure  2.  Probability  histogram  (upper)  and  autocorrelation  function  (lower)  of  the  surface 
height  fluctuation  as  calculated  from  Talystep  measurements  for  the  surface  ( #46)  used  in 
this  paper. 
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light  was  measured  at  four  wavelengths  (0.63  pm,  1  15  /an,  3.39  /an  and  10.6  /an)  using 
the  equipment  described  in  [10].  For  each  angle  of  incidence,  measurements  are  made 
with  p-polarization  incident  and  p-polarization  collected  (‘p-p’  scattering)  and  s- 
polarization  incident  and  s-polarization  collected  (‘s-s’  scattering);  no  crossed  polarized 
components  were  detectable.  For  a  perfect  conductor,  these  measurements  give  a 
complete  description  of  the  scattering  characteristics  of  the  surface,  but  in  general  four 
scattering  coefficients  are  required  for  materials  of  finite  conductivity;  these  can  be  found 
by  measuring  the  polarization  of  the  scattered  light  for  various  input  polarizations.  Also, 
the  measurements  reported  here  yield  the  relative  scattering  cross  section,  as  no  absolute 
calibration  is  performed. 

The  relative  scattering  cross  sections  for  angles  of  incidence  of  0°,  -30°  and  —60° 
and  the  four  wavelengths  are  shown  in  figures  3  and  4  for  the  gold-coated  surface  and 
figures  5  and  6  for  the  dielectric  surface.  The  enhanced  backscatter  peak,  where  present, 
occurs  on  the  right-hand  side  of  the  graphs  (i.e.  at  positive  angles,  see  figure  1  for  the 
sign  convention  for  the  angles)  and  any  specular  component  is  on  the  left-hand  side 
(i.e.  negative  angles);  for  the  10.6  /im  measurements,  the  specular  peak  was  very  much 
greater  than  the  diffuse  component  and  is  not  shown.  A  few  features  are  of  particular  note. 
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Figure  3.  Relative  scattering  cross  section  as  a  function  of  scattering  angle  for  the  gold-coated 

surface,  for  angles  of  incidence  of0°,  —30°  and  —60°,  for  p-p  scattering  (open  circles)  and  i 

s-s  scattering  (crosses).  The  left-hand  column  is  for  a  wavelength  of  0.63  /an,  for  which 

a/X  *  1.93  and  z/X  -  5.02  and  the  right-hand  column  for  2  «=  1.15  /tm,  for  which  a/X  =  1.07 

and  z/X  =■  2.76.  The  enhanced  backscatter  peak,  where  present,  occurs  at  positive  angles 

(right-hand  side  of  each  graph). 
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(1)  The  enhanced  backscatter  peak  and  sidelobe  structure  are  clearly  visible  for  the 
shorter  wavelengths  at  an  angle  of  incidence  less  than  approximately  —30°  for  the 
gold-coated  surface;  the  width  of  the  peak  is  proportional  to  the  wavelength.  The  peak 
is  not  observed  for  the  p-p  scattering  at  10.6  /an  for  the  gold  surface  for  scattering 
from  the  dielectric. 

(2)  With  the  exception  of  the  p-p  case  at  10.6  /an,  the  p-p  and  s-s  scattering  by  the 
gold  surface  are  very  similar;  for  the  dielectric  surface,  however,  the  p-p  and  s-s  scattering 
cross  sections  are  quite  different,  as  one  might  expect  by  analogy  with  reflection  from 
a  planar  surface.  Using  a  value  of  n  =  1.41  for  the  refractive  index  of  the  (almost  lossless) 
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Figure  4.  As  for  figure  3  but  wavelengths  of  3.39  fim  ( a/X  =  0.36,  x/X  =  0.94)  and  10.6 pm 
(a/X  «=  0.12,  x/X  =  0.30).  For  the  10.6  pan  curves,  the  angles  of  incidence  were  0°,  -20°  and 
—40°..  The  (strong)  specular  component  in  the  10.6  pan  curves  is  not  shown 


dielectric  gives  a  Brewster  angle  of  «  55  °.  Considering  single  scattering  to  be  the  dominant 
mechanism  and  treating  this  as  a  reflection  from  a  locally  plane  surface  gives  an 
expected  minimum  of  the  p-p  scattered  intensity  at  an  angle  equal  to  approximately 
( - 110°  -  incident  angle):  the  angles  are  roughly  in  accordance  with  this  simple  picture. 
The  s-s  and  p-p  scattered  intensities  in  the  backscatter  direction  appear  to  be  almost 
equal  to  each  other  for  all  angles  of  incidence  and  wavelengths,  for  the  dielectric. 

(3)  The  overall  shape  of  the  curves  is  dramatically  different  from  the  Gaussian-type 
shapes  (centred  on  the  specular  angle)  normally  encountered  in  scattering  from 
low-sloped  surfaces. 

The  principal  purpose  of  figures  3  to  6  is  to  provide  a  reliable  set  of  data  for 
comparison  with  numerical  calculations,  and  analytical  theories  should  any  become 
available. 
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Figure  5.  As  for  figure  3  but  for  the  dielectric  surface,  n0  43  -  -  1.41,  n,  15  =  1.40.  The  vertical 
scale  is  not  the  same  as  that  used  in  figure  3  (both  are  relative  scattering  cross  sections). 


3.  Kirchboff  multiple  scatter  approximation 

One  can  compare  the  above  experimental  results  to  those  of  ‘  exact’  numerical  calculations 
based  on  the  extinction  theorem  and  its  extensions  [3-7],  and  some  comparisons  of 
experiment  and  calculations  for  a  perfect  conductor  were  given  in  [  10].  Although  such 
comparisons  are  valuable,  one  problem  with  the  ‘exact’  numerical  solution  is  that  it 
gives  little  physical  insight  into  the  problem.  We  have  therefore  attempted  to  extend 
the  KirchhofF approximation  (i.e.  tangent  plane  approximation  for  each  scattering  event) 
to  double  (and  multiple)  scattering  [11]. 

The  numerical  calculations  were  carried  out  using  the  method  described  in  [11]  for 
a  perfect  conductor;  typically  the  energy  conservation  (unitarity)  held  to  better  than 
3%  considering  just  the  single  and  double  scatter  terms  for  surface  #46  (except  for  the 
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Figure  6.  As  for  figure  4  but  for  the  dielectric  surface,  n33,  =  nl0  6  =  1.41. 


— 60°  results  for  which  the  departure  from  unitarity  is  «  6  % )  and  the  results  are  averaged 
over  approximately  103  realizations.  Figures  7  and  8  show  the  results  of  the  calculations 
for  A  —  0.63  ftm  and  1.15  /tin  respectively,  for  incident  angles  of  0°,  —  30°  and  —60°  and 
s-s  aud  p-p  scattering.  Each  graph  shows  the  single,  double  and  total  scattered  intensity. 
The  enhanced  backscatter  peak  occurs  only  in  the  double  scatter  component,  showing 
conclusively  that  the  enhancement  is  a  multiple  scattering  effect  The  enhancement  is 
on  the  order  of  a  factor  of  two  in  the  double  scattered  component  for  all  angles  of 
incidence,  but  the  enhancement  in  the  total  intensity  is  much  less  than  two  and  decreases 
with  increasing  angle  of  incidence  due  to  the  fact  that  the  double  scattered  intensity 
also  decreases  with  incidence  angle 

Figure  9  compares  the  total  scattered  intensity  for  s-s  scattering  from  figures  7  and 
8  with  the  results  of  ‘exact’  numerical  calculations  (based  on  the  extinction  theorem 
method  for  a  perfect  conductor  [3])  and  the  experimental  results  of  figure  3,  for 
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s-s 


p-p 


Figure  7.  Numerically  calculated  scattered  intensities  for  a  perfect  conductor,  using  the 
double-scattering  Kirdihoff  approximation,  for  angles  of  incidence  of  0°,  —30°  and  —60° 
and  s-s  scattering  (left)  and  p-p  scattering  'right).  The  wavelength  is  0.63 /an.  Each  graph 
shows  the  doubly  scattered  intensity  (lowest  curve),  single  scattered  (middle)  and  total 
intensity  (coherent  sum)  (top  curve). 


X  —  0.63  job  and  1.15  pm  at  three  angles  of  incidence.  The  two  numerical  calculations 
agree  well,  showing  that  the  Kirchhoff  approximation  is  reasonable  for  these  surface 
parameters  (the  average  radius  of  curvature,  defined  as  the  inverse  of  the  standard 
deviation  of  the  surface  curvature  2n/3<t/t2,  is  «2.4/tm  for  surface  #46)  and  both 
agree  well  with  the  experimental  measurements  for  zero  angle  of  incidence.  However, 
there  is  a  clear  discrepancy  between  experiment  and  numerical  calculation  for  the  -30° 
and  —60°  angles  of  incidence.  (This  general  behaviour  is  also  shown  in  the  case  of  p-p 
scattering.) 
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Figure  8.  As  for  figure  7,  but  wavelength  of  1.15  /on. 
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One  possible  cause  of  the  discrepancy  could  be  that  the  calculations  are  for  a  perfect 
conductor,  whereas  the  experiments  are  for  a  real  metal  (gold).  However,  calculations 
by  ourselves  and  others  [5]  show  that,  for  these  values  of  rms  surface  height  and 
correlation  length,  there  is  very  little  difference  between  the  results  for  gold  and  for  a 
perfect  conductor,  particularly  for  s-s  scattering.  One  problem  with  most  methods  of 
calculation,  including  that  used  here,  is  that  a  very  small  length  of  surface  is  vOncicLred, 
giving  rise  to  the  possibility  of  an  ‘end-effect’  error  (e.g.  due  to  long-range  surface 
plasmons);  however,  the  method  of  calculation  of  Saillard  and  Maystre  [7]  uses  an 
extremely  long  length  of  surface  with  good  agreement  with  the  other  calculations  and 
poor  agreement  with  the  measurements  at  larger  angles  of  incidence. 

It  seems,  therefore,  that  there  may  be  some  aspect  of  the  experiment  that  does  not 
correspond  to  the  calculations.  Previous  results  for  a  Lambert>an  diffuser  have 
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Figure  9.  Comparison  of  Kirchhoflf  calculation  (from  figures  7  and  8),  ‘exact’  numerical 
calculation  and  experiment  (from  figure  3),  for  s-s  scattering  at  3  =  0.63  fan  and  1.15  pm, 
and  angles  of  incidence  equal  to  O',  —30°  and  —60°.  The  solid  curves  are  the  Kiichhoff 
calculation,  triangles  the  ‘exact’  calculation  and  crosses  are  the  experimental  results.  Note 
the  good  agreement  between  the  two  numerical  calculations  but  the  departure  of  the 
experimental  results  for  larger  angles  of  incidence. 


demonstrated  that  the  scatterometer  measures  the  correct  quantity  [10].  The  measure¬ 
ment  of  surface  properties  might  be  in  error.  If  one  calculates  the  scattered  intensity 
for,  say,  -60°  angle  of  incidence  for  a  surface  that  has  an  rms  roughness  50%  larger 
than  the  measured  value,  then  reasonable  agreement  is  obtained  between  experiment 
and  numerical  calculation.  However,  (a)  it  is  extremely  unlikely  that  such  a  gross  error 
could  occur  (stylus  tips  effects  are  discussed  by  Church  [12])  and  (b)  the  agreement 
for  0°  angle  of  incidence  is  then  very  poor  indeed,  particularly  as  regards  the  location 
of  the  minima  around  the  backscatter  peak.  Ishhnaru  and  Chen  [13]  have  shown  that 
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a  departure  from  Gaussianity  of  the  correlation  function  could  be  responsible  for  the 
discrepancy,  and  the  measured  coi  relation  does  show  a  small  departure  from  the 
Gaussian  shape.  However,  it  is  notoriously  difficult  to  estimate  the  correlation  function 
of  stylus  traces  and  the  departure  shown  in  figure  2  is  characteristic  of  inadequate 
de-trending  of  the  mean;  the  method  of  manufacture  of  the  surfaces  strongly  encourages 
a  Gaussian  correlation  of  surface  height.  The  cause  of  this  discrepancy  for  larger  angles 
of  incidence  is  therefore  not  resolved  at  tire  present  time. 


4.  Summary 

A  set  of  scattering  data  for  a  one-dimensional  surface  at  four  wavelengths,  three  angles 
of  incidence  and  two  materials  has  been  presented  for  critical  comparison  with  numerical 
calculations  and  theoretical  studies.  A  multiple  scatter  extension  of  the  Kirchhoff 
approximation  has  been  shown  to  provide  additional  physical  evidence  that  the 
predominant  cause  of  the  enhanced  backscatter  peak  is  due  to  multiple  scattering.  There 
remains  a  significant  disagreement  between  experiment  and  numerical  calculations  for 
large  angles  of  incidence  the  cause  of  which  is  still  unresolved. 

The  data  presented  in  figures  3  to  6,  together  with  sample  Talystep  traces,  is  available 
on  a  PC-  or  Macintosh-compatible  diskette  on  application  to  the  first  author. 
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The  Kirchhoff  approximation  is  used  to  find  the  scattered  intensities  and  the  Mueller  matrix  for  a  gaussian  surface  with  stan¬ 
dard  deviation  of  height  a— 1.22  pm  and  correlation  length  t=3.I8  pm  using  a  wavelength  of  2=0.633  pm.  The  cases  of  gold 
(n=0.167+i  3.149)  and  dielectric  (n=  1.41)  surfaces  are  considered.  Results  are  compared  with  experimental  observations  and 
are  shown  to  give  good  agreement.  The  single  and  double  scattered  intensities  are  shown  to  be  fully  polarised  for  both  cases 
whereas  the  total  (single  plus  double)  shows  a  strong  unpolarised  component  for  the  gold  case.  This  effect  is  understood  by 
considering  that  the  single  and  double  scatter  terms  ate  partially  coherent  and  at  right  angles  to  each  other. 


1.  Introduction 

Scattering  of  electromagnetic  waves  from  one-di¬ 
mensional  rough  surfaces  has  been  a  topic  of  much 
interest  in  recent  years,  the  motivation  for  a  great 
deal  of  this  work  being  the  phenomenon  of  enhanced 
backscatter  [1,2].  Experimental  [3-5]  and  theoret¬ 
ical  [6-10]  curves  for  die  distribution  of  scattered 
intensity  for  the  two  cases  of  s-polarisation  (TE)  in¬ 
cident  and  scattered  (denoted  by  ss)  and  p-polar- 
isatidn  <TM)  incident  and  scattered  (denoted  by  pp) 
have'been  found.  These  cases  have  been  considered, 
since  for  a  one-dimensional  surface  ( i.e.  a  surface  for 
which’tlre  height  is  constant  in  one  direction)  there 
is  no  scattering  from  s  to  p  or  from  p  to  s;  s  and  p 
are  the  two  fundamental  polarisations  into  which  all 
other  polarised  cases  can  be  decomposed.  However, 
as  has  been  pointed  out  recently  [5],  the  ss  and  pp 
intensity  distributions  do  not  give  a  fall  represen¬ 
tation  of  die  scattering  process  for  these  surfaces.  This 
is  due  to  the  fact  that  the  relative  phase  of  the  ss  and 
pp  scattered  fields  is  not  contained  in  the  ss  and  pp 
;  intensities.  For  example  an  incident  field  polarised 
■  at  45°  to  the  plane  of  incidence  may  give  a  com- 
{  ponent  of  the  scattered  field  circularly  polarised  (i.e. 
«  an  incident  field  with  s  and  p  components  the  same 
|  amplitude  and  in  phase  could  give  a  scattered  field 


with  s  and  p  the  same  amplitude  with  a  n/2  phase 
difference).  In  the  ss  and  pp  intensities  these  phase 
terms  are  not  used  since  the  modulus  squared  of  the 
scattered  field  is  calculated. 

In  calculations  of  the  intensity  distribution  the 
scattered  field  is  first  evaluated  and  the  modulus 
square  then  taken.  This  means  that  the  phase  term 
can  be  found  directly  from  the  calculations.  How¬ 
ever,  in  experiments  only  the  intensity  can  be  mea¬ 
sured  so  the  phase  difference  cannot  be  directly 
found.  Instead  the  scattered  intensities  for  certain 
combinations  of  incident  and  detected  polarisations 
are  required  to  find  the  Stokes  parameters  of  the 
scattered  light  [11].. 

In  this  communication  the  recently  developed 
multiple  scattering  Kirchhoff  approximation  (KA) 
[  12-16]  is  used  to  calculate  the  Stokes  parameters 
for  a  particular  surface.  The  advantage  of  the  KA  is 
that  the  single  and  double  scatter  terms  can  be  sep¬ 
arated  i.e.  it  is  possible  to  separate  light  paths  inter¬ 
secting  the  surface  profile  once  or  twice.  The  cal¬ 
culation  has  been  performed  for  a  gold  surface  with 
refractive  index  n=0.167+i3.149  (at  2=0.633  pm) 
and  a  dielectric  surface  of  n= 1.41.. 

The  theory  of  the  calculations  and  experimental 
procedure  is  discussed  in  sect.  2  and  comparisons  of 
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numerical  and  experimental  results  are  given  in  sect. 


'i 


3. 


2.  Theory 

In  the  Kirchhoff approximation  the  total  field  on 
the  surface  is  approximated  by  the  incident  field  plus 
the  reflected  field  at  each  point.  The  reflection  is  as¬ 
sumed  to  be  from  a  plane  to  permit  Fresnel  reflec¬ 
tion  coefficients  to  be  used.  Tins  gives  a  condition 
for  the  approximation  to  be  valid  -  the  radius  of  cur¬ 
vature  at  each  point  of  the  surface  must  be  large 
compared  to  the  illuminating  wavelength.  The  single 
scatter  term  is  given  by  [  1 2, 1 3  ] 

./ 2  V/2<aP (ip)  l+cosffl+ft) 

*  '  \xkr)  2  costf+cos#, 

xjsfx'.z'^explif^-iiJ-^ldjr',  (1) 

r 

where  6i  and  0  are  the  incident  and  scatter  angles  re¬ 
spectively,  <f  is  a  phase  term  independent  of  x\  R  is 
the  position  vector  of  the  point  (x\  z’),  R  is  the 
planar  reflectivity  for  the  tangent  at  (x1,  z')  and  S(z', 
z')  is  a  shadow  function  depending  on  both  0,  and 
0.  The  double  scatter  contribution  can  be  considered 
as  two  terms,  first  the  field  at  any  surface  point  due 
to  light  from  one  other  surface  point  is 

E,(xi, z2)  =  -  S.*|(  1  +*i )[*",  ) 

ra  J 

+cos0,)tfj1)(fa-12)j£1(.ir1,z1) ,  (2) 

with  Sa  a  shadow  function  depending  on  6,  and  the 
vector  between  points  1  and  2,  m,  the  gradient  of  the 
surface  at  the  first  point,  Rt  the  Fresnel  reflection 
coefficient  at  the  first  point  and  Hj>l}(krl2)  and 
//f  ’  )(/cr,2 )  the  zeroth  and  first  order  Hankel  func¬ 
tions  of  the  first  kind  The  total  field  due  to  all  com¬ 
binations  of  first  and  second  points  is 


(j&)'  S*4ti 

v  r  r 


Xfc(m2sin0-cos0)£  (x2,z2) 

+  i(l-.ft2)  |  exp(iA-f?2)dxldx2 ,  ' 

07*2  J 

(3)  , 

where  S'  represents  the  shadow  function  for  the  scat-  j 

tered  radiation,  m2  is  the  gradient  of  the  surface  at 
the  second  point  and  R2  is  the  Fresnel  reflection 
coefficient  at  th-  second  surface  point.  The  shadow 
functions  used  here  are  geometric  (straight  line) 
shadow  functions  and  were  calculated  explicitly  for 
every  point  (or  combination  of  prints)  by  perform¬ 
ing  a  ray  trace.  If  the  ray  of  interest  was  intersected 
by  the  surface  at  any  point  the  shadow  function  took 
the  value  zero  and  the  contribution  'rom  that  point 
was  neglected.  Eqs.  (2)  and  (3)  tog  ther  represent 
the  double  scatter  contribution.  Any  material  de¬ 
pendence  is  entirely  contained  within  toe  Fresnel  re¬ 
flectivity  coefficients  R,  Rl  and  R2.  Ligh'  paths  trav¬ 
ersing  inside  ihe  material  are  not  induled  in  this 
method  but  it  is  expected  that  the  contribution  from 
such  paths  is  small  since  the  necessary  scatter  angles 
are  large  compared  to  the  refraction  angles  of  the 
materials  considered.  The  simplest  case  occurs  for  a 
perfect  conductor  for  which  R- 1  for  p-polarisation 
and  R=  - 1  for  s  causing  most  of  the  above  teitnsto 
drop  out,  leaving  only  one  contribution  for  each  case. 

Eqs.  (1),  (2)  and  (3)  are  discretiscd  and  oro- 
grammed  into  a  computer.  Intensity  curves  are  .cal¬ 
culated  for  particular  surface  profiles  and  the  inten¬ 
sities  averaged  over  typically  800  different  surfaces  \ 

of  the  same  statistics  to  produce  the  final  curves. 

The  convergence  of  the  series  given  by  the  sum  of 
the  Kirchhoff  single  scatter  plus  double  scatter  plus  I 

etc.  has  not  been  proved  rigorously.  However,  cal¬ 
culations  of  the  triple  scatter  term  show  that  the  in¬ 
tegrated  energy  in  this  order  of  scattering  is  much 
smaller  than  in  the  double  scatter  (approximately 
1  /50th  the  double  scatter).  While  this  does  not  prove 
the  convergence  of  the  series,  it  at  least  shows  a  re¬ 
quired  effect  -  rhe  relative  importance  of  the  results 
of  higher  order  calculations  should  decrease  as  the 
order  increases. 

The  scattering  of  the  field  can  be  denoted  in  ma-  1 

trix  terms  as 


i 

- 1 
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where  E#  and  £,  are  the  scattered  p  and  s  fields,  £ip 
and  E*  are  the  incident  p  and  s  fields  and  and  .4 
are  the  field  scattering  coefficients.  This  can  be  writ¬ 
ten  since  any  polarised  state  can  be  decomposed  into 
s  and  p  components  and  there  is  no  depolarisation 
for  scattering  of  s  and  p  from  one-dimensional  sur¬ 
faces.  Note  that  all  the  terms  represent  complex 
numbers  and  are  functions  of  angles.  Hie  Stokes  vec¬ 
tor  is  given  by  four  parameters  which  give  for  the 
scattered  wave  (11) 

/=<£p£;  +£,£?>,  <£„£;—£,£?>  , 
£/=<£„£?+£,£;>,  T=i<£p£j  —  £,£p>  , 

(5) 

where  <  >  denotes  the  average  over  different  sur¬ 
faces  of  the  same  statistics.  These  give  a  full  char¬ 
acterisation  of  the  polarisation  properties  of  the  scat¬ 
tered  wave. 

The  Stokes  parameters  of  the  scattered  field  are 
related  to  those  of  the  incident  field  by  the  Mueller 
matrix,  M 


S—MSo , 


(6) 


where  Mis  a  4x4  matrix  which,  for  the  one-dimen¬ 
sional  case  is  given  by 


(mH  mI2  0  0  ' 

ml2  mu  0  0 

0  0  mj* 

0  0  —  m34  si 


(7) 


This  matrix  describes  fully  the  scattering  properties 
of  the  target  surface.  For  an  incident  polarisation  at 
+45°  (halfway  between  the  s  and  p  polarisation  di¬ 
rections)  the  incident  Stokes  vector  is 


so  that  the  final  Stokes  vector  of  the  scattered  light 
is,  from  e*;.  (6) 


/  \ 


S=\ 


mss 


.  \-mMJ 


(9) 


i.e.  the  final  Stokes  vector  is  just  given  by  the  sep¬ 
arate  non-zero  Mueller  matrix  elements.  In  the  cal¬ 
culations  £ip  and  £„  are  set  to  1/-J2  to  give  an  in¬ 
cident  unit  amplitude  plane  wave  linearly  polarised 
at  +45".  The  scattered  s  and  p  fields  are  found  from 
eqs.  ( 1 )  to  (3 )  and  the  Stokes  matrix,  and  hence  the 
Mueller  elements,  are  found  by  using  eq.  (5). 

Experimentally  the  Mueller  matrix  elements  are 
found  by  measuring  the  scattered  intensity  distri¬ 
bution  using  suitable  combinations  of  incident  and 
detected  polarisations  [5,11,14].  The  terms  m , ,  and 
m,2  are  measured  by  using  p-poiarisation  incident 
and  p  detected  (given  by  /„)  and  s  incident  and  s 
detected  (/,).  The  signals  resulting  are  added  for  m, , 
and  subtracted  (p-polarised-s-polarised)  for  ml2. 
mss  is  found  by  scattering  light  polarised  at  +45° 
then  finding  the  difference  of  the  scattered  light  po- 
larisedat  +45°  (/+)  minus  that  polarised  at  -45° 
(/_ ).  The  final  term  is  found  by  again  using  incident 
light  polarised  at  +45°  and  finding  the  difference  of 
left  circularly  scattered  light  (/L)  minus  right  cir¬ 
cularly  scattered  light  (/R)  [17].. 

The  four  Stokes  parameters  obey  the  following  re¬ 
lation  [18,19] 

I2ZQ2+U2+V2,  (10) 

the  equality  holding  for  fully  polarised  light.  For  the 
case  of  incident  polarisation  at  +45°  the  terms  of 
this  expression  are  given  by  eq.  (9).  It  is  thus  pos¬ 
sible  to  separate  the  polarised  and  unpolarised  com¬ 
ponents  of  the  scattered  energy  for  this  situation,  the 
polarised  being  given  by  (Q2+U2+  V2),/2  and  the 
unpolarised  by  /—  (Q2  +  C/2+V2),/2. 


3.  Results 

The  validity  of  using  the  Kirchhoff  approximation 
is  limited  to  surfaces  containing  radii  of  curvature 
greater  than  the  illuminating  wavelength  to  permit 
the  use  of  planar  reflection  coefficients  in  eqs.  ( 1 ), 
(2)  and  (3).  The  applicability  of  the  method  to  the 
surface  used  can  be  verified  by  noting  that  the  radius 
of  curvature  is  approximately  the  inverse  of  the  sec¬ 
ond  derivative  of  the  surface  profile.  As  the  height 
distribution  is  gaussian  so  is  the  second  derivative 
with  a  standard  deviation  ah~  given  by 
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where  oj,  is  the  standard  deviation  of  the  gaussian 
probability  distribution  of  heights  and  t  is  the  1/e 
level  of  the  gaussian  autocorrelation  function  for  the 
surface  of  interest  This  gives  a  value  of  <rh.s=0,42 
urn- 1  for  the  case  considered  here  (surface  #46  with 
a- 1.22  pm  and  r=3.18  pm).  Therefore,  approxi¬ 
mately  67%  of  the  radii  of  curvature  will  have  values 
greater  than  2.4  pm,  so  most  of  the  curvatures  have 
values  greater  than  the  wavelength  leading  to  the 
conclusion  that  it  is  reasonable  to  expect  the  Kir- 
chhoff  approximation  to  be  valid  in  this  case. 

The  calculated  unitarity  or  relative  scattered  en¬ 
ergy  for  the  gold  and  dielectric  surfaces  are  given  in 
table  1  for  the  single,  double  and  total  (single  plus 
double  added  coherently)  components  of  the  fun¬ 
damental  polarisations  s  and  p.  The  values  of  total 
scattered  energy  for  these  surfaces  are  less  than  one 
due  to  the  finite  conductivity  of  the  surface  materials 
causing  light  to  be  transmitted  or  absorbed. 

Figs.  1  and  2  show  the  Kirchhoff  single  and  double 
scatter  terms  for  the  gold  and  dielectric  cases  re¬ 
spectively.  Note  that  there  are  different  scales  on  the 
single  and  double  graphs  for  both  situations.  The 
vertical  dashed  line  on  the  graphs  marks  the  incident 
direction  so  that  light  scattered  in  this  direction  is 
backscattered.  Enhanced  backscatter  effects  can  be 
seen  in  the  double  scatter  terms  only,  as  expected 
from  previous  work  [12-16].,  In  the  gold  case  the 


Table  1 
Gold  surface 


ml2  term  is  very  small  for  both  single  and  double 
which  means  that  there  is  very  little  difference  be¬ 
tween  the  s  and  p  scattered  intensities.  In  the  die¬ 
lectric  curves  this  difference  is  more  important,  as 
expected,  since  the  Fresnel  reflection  coefficients 
show  a  bigger  difference  for  a  dielectric  than  for  a 
good  conductor. 

The  total  (single  and  double)  scattered  intensities 
are  shown  in  figs.  3  and  4  for  the  gold  and  dielectric 
cases  respectively  and  are  compared  with  the  exper¬ 
imentally  measured  values.  In  the  gold  case  the 
agreement  is  reasonably  good  at  low  angles  of  inci¬ 
dence  apart  from  the  term  which  is  much  smaller 
in  the  calculations.  As  in  other  measurements  [8,12] 
the  agreement  at  higher  angles  of  incidence  is  not  so 
good;  this  appears  to  be  the  case  for  all  methods  of 
calculation,  not  just  the  Kirchhoff  method.  The  rea¬ 
son  for  the  difference  in  the  ms<  term,  which  is  the 
difference  between  right  and  left  circularly  polarised 
light,  may  be  that  the  refractive  index  of  the  material 
is  not  exactly  as  used  in  the  calculation  so  that  the 
phase  difference  between  the  s  and  p  scattered  com¬ 
ponents  is  different  to  that  expected.  The  agreement 
for  the  dielectric  case  is  much  better.  Since  the  re» 
fractive  index  is  real  for  this  material  there  is  no 
scattering  from  linear  to  circular  polarisation  so  m34 
is  near  zero  in  this  case. 

The  degree  of  polarisation  of  the  scattered  light  is 
illustrated  in  figs.  5  and  6  for  the  calculations  of  the 
single  and  double  sc  *ter  terms  and  figs.  7  and  8  for 
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Fig-  5.  Graphs  of  the  polarised  and  unpolarised  scatter  components  as  a  function  of  scatter  angle  for  0“  incidence  (top),  20°  incidence 
(middle)  and  40°  incidence  (bottom).  Single  scatter  (left )  and  double  scatter  (right )  for  gold  coated  surface  #46. 
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the  total  and  experimental  curve*.  In  both  the  gold 
and  dielectric  cases  thy  single  and  double  scattered 
intensities  (figs.  5  and  6)  are  fully  polarised  whereas 
in  the  total  and  experimental  curves  (figs.  7  and  8) 
the  gold  shows  a  large  unpolarised  component  and 
the  dielectric  case  is  again  almost  fully  polarised.  It 
is  believed  that  this  ts  caused  by  the  partial  coher¬ 
ence  of  the  single  and  double  scattered  fields.  Figs. 
1  and  2  show  that  m33  is  negative  for  single  scatter 
and  positive  for  double  scatter  for  both  materials. 
The  experimental  interpretation  of  m33  is  that  it  rep¬ 
resents  the  difference  in  the  intensity  polarised  at 
t- 45°  and  the  intensity  at  -45°  for  incident  polar¬ 
isation  at  +45  This  means  that  the  single  and  dou¬ 
ble  scattered  fields  have  linearly  polarised  compo¬ 
nents  at  right  angles  to  one  another.  For  scattering 
from  a  perfect  conductor  the  singly  scattered  field  is 
the  same  polarisation  &>  the  incident  field  but  the 
doubly  scattered  is  orthogonally  polarised  to  the  in¬ 
cident  (fig.  9).  Also  note  that  the  total  curves  can  be 
obtained  simply  by  summing  the  single  and  double 
curves  i.e.  the  single  and  double  terms  are  largely  in¬ 
coherent  f  12—16].  Two  incoherent  components  at 
right  angles  give  a  resultant  field  which  has  a  varying 
polarisation  (the  phase  between  the  two  components 
is  not  fixed )  f  18]  and  so  an  unpolarised  component 
results.  It  is  important  to  note  that  the  single  and 
double  scattered  intensities  are  fully  polarised  but 
because  they  have  components  polarised  in  different 
directions  and  are  incoherent  they  produce  unpolar¬ 
ised  light  when  combined.  In  the  calculations  the  un- 
polansed  component  contains  the  backscatter  peak 
since  the  enhancement  is  contained  in  the  double 
scattered  component  and  so  gives  an  unpolarised  re¬ 
sult  when  combined  with  the  single  scatter  term.  The 
same  effect  will  occur  in  scattering  from  two-dimen¬ 
sional  surfaces  for  scattering  in  the  plane  of  inci¬ 
dence  although  a  simple  geometrical  effect  could  ac¬ 
count  for  the  unpolarised  component  out  of  this 
plane.  This  effect  is  not  seen  in  the  dielectric  curves 
simply  because  the  double  scatter  is  a  much  smaller 
effect  for  this  case.  In  table  1  the  double  scatter  en¬ 
ergy  is  always  one  or  two  orders  of  magnitude  down 
on  the  single  scatter  for  the  dielectric  surface  so  this 
polarisation  effect  will  not  be  visible. 

In  fig.  7  the  experiment?.!  curve  shows  a  peak  in 
the  polarised  component  which  is  not  in  the  calcu¬ 
lated  values.  This  may  be  due  to  the  effect  of  triple 
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Fig.  9.  Figure  showing  +45*  incident  polarisation  and  resultant 
single  and  double  scatter  polarisations  lor  a  perfect  conductor. 
The  dots  represent  an  electric  field  out  of  the  plane  of  the  dia¬ 
gram,  the  crosses  an  electric  field  into  the  plane.  In  the  boxes  at 
the  top  of  rite  diagram  an  arrow  downwards  represents  an  s  field 
out  of  the  plane.  Note  the  single  scatter  is  co-polarised  with  the 
incident  and  the  double  scatter  is  cross-polanscd, 

scatter  which  gives  a  term  polari'ied  in  the  same  di- 
rectic.i  as  the  single  scatter.  Even  if  these  terms  are 
incoherent  the  plane  of  polarisation  will  not  vary. 
Therefore,  since  the  triple  scatter  is  expected  to  show 
an  enhanced  backscatter  peak,  a  peak  in  the  pola¬ 
rised  component  will  result.  The  triple  scatter  can  be 
expected  to  decrease  as  the  angle  of  incidence  is  in¬ 
creased  (cf.  the  double  scatter  term  in  refs.  [12]  and 
[13])  so  that  the  experiment  and  theory  become 
more  compara'  !e  for  higher  angles  of  incidence. 


4.  Conclusions 

The  Mueller  matrix  for  a  particular  surface  has 
been  calculated  using  the  Kirchhoff  method  which 
separates  the  single  and  double  scatter  contributions. 
The  advantage  of  using  this  method  is  seen  when  the 
difference  between  the  polarised  and  unpolarised 
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components  are  calculated.  The  single  and  double 
scatter  terms  are  fully  polarised  but  the  sum  of  <he 
two  is  partially  polarised  for  a  gold  coated  surface 
since  it  is  made  >ip  of  the  addition  of  two  partially 
coherent  components  at  right  angles.  This  effect  is 
not  seen  in  a  dielectric  surface  because  the  double 
scatter  is  too  weak. 
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THEORY  AND  MEASUREMENTS  OF  THE  ANGULAR 
CORRELATION  OF  SPECKLE  PATTERNS 

R  W  Syratt,  J  C  Dainty 
Imperial  College,  London 

In  this  work,  we  illuminate  a  randomly  rough  surface  at  one  angle,  6\ ,  record 
the  intensity  around  a  small  range  of  angles  centered  on  another  angle  9\  and 
repeat  this  process  for  another  angle  of  illumination  and  viewing,  9  and  9' . 
The  two  recorded  speckle  patterns  are  then  cross-correlated  to  obtain  the 
correlation  coefficient  between  them. 
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Figure  1 
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Figure  2 


The  work  conducted  so  far  is  limited  to  observation  in  the  specular  direction 
( 9 ;  =  9j,  6’  =  9)  this  ensures  a  large  correlation  range,  but  limits  the 
technique  to  examination  of  surfaces  which  do  not  reflect  a  large  specular 
component.  The  geometry  can  be  seen  in  Figure  1.  A  range  of  speckle 
patterns  (changing  59)  is  correlated  with  a  fixed  reference  pattern  (9i). 

Figure  2  schematically  shows  the  experimental  equipment;  typically  twenty- 
five  speckle  images  are  taken  from  the  charge-coupled-device  (CCD)  camera 
and  averaged  to  reduce  measurement  noise.  The  optical  disk  is  necessary  to 
store  the  vast  amounts  of  data  gathered. 
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Figure  3 

Graphs  of  experimental  angular  correlation  curves  at  different  references 
angles,  for  a  Gaussian  random  rough  surface  —  with  an  rms  height  (a)  of 
2. 27 ± 0.2pm  and  correlation  length  of  20. 7+ 0.2pm  —  illuminated  by  a  HeNe 
laser  of  wavelength  (A)  0.633  pm  are  given  in  Figure  3;  these  are 
computationally  calculated  from  the  speckle  images.  Overlaying  them  are  the 
curves  from  an  equation  derived  from  the  single  scatter  Kirchoff 
approximation,  a  form  of  which  is  given  by  equation  1.  In  practice  a  more 
general  two  dimensional  equation  is  used  and  averaged  over  the  field  view  of 
the  CCD  camera, 

c(e-  “)=5^Tsij  “p  (“se‘'  ros!(f,'+M»2  ]  (D 
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SCATTERING  FROM  DENSE  VOLUMES 


N  C  Bruce 

Imperial  College,  London 
Introduction 

There  has  recently  been  a  great  deal  of  interest  in  the  use  of  light  for  non- 
invasive  testing  in  medicine.  Possible  applications  include  mammography, 
investigation  of  skin  cancers  and  detection  of  cancerous  tissue  during  surgery. 
The  use  of  optics  in  a  clinical  environment  requires  compact,  robust  devices 
which  are  possible  by  using  fibre  optics.  An  experimental  model  of  such  a, 
device  has  been  built  and  its  response  for  different  experimental  conditions 
studied. 


Experiment 

The  experimental  arrangement  is  shown 
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Light  from  a  He-Ne  laser  is  delivered  to  the  sample  via  a  mono-mode  fibre  to 
give  a  clean  illumination  spot  in  the  sample.  A  multi-mode  fibre  is  used  to 
collect  the  scattered  light  and  deliver  it  to  a  photon  counting  photo-multiplier 
connected  to  a  computer.  The  detection  fibre  is  scanned  away  from  the 
source  fibre,  always  in  the  same  plane.  The  sample  holder  was  a  stainless  steel 
pot  of  2 cm  diameter.  This  pot  was  large  enough  so  that  the  effect  of 
reflections  from  the  edge  was  negligible.  The  samples  used  were  collections 
of  polystyrene  spheres  in  water  obtained  from  Sigma  Chemical  Company. 
The  samples  were  10%  by  volume  concentration  of  particles  of  diameters 
0.091/im,  0.296 fim,  0.4 6pim„  0.605 jim  and  0.11Zjj.m.  Water  was  added  to 
some  of  these  samples  to  give  1%  volume  density  liquids.  Since  these  samples 
are  so  well  characterised  it  will  be  possible  to  compare  the  experimental 
results  with  theoretical  values. 


Results 

The  figures  below  show  the  number  of  counts  per  second  versus  separation  of 
the  two  fibres  for  the  10%  and  1%  volume  density  samples  for  the  0.091  jim. 
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From  these  curves  it  can  be  seen  that  there  is  no  linear  relationship  between 
the  width  or  the  maximum  value  of  the  scatter  pattern  and  the  size  of  the 
particle.  Indeed  from  the  diffusion  theory  [1]  the  figure  that  describes  the 
scattering  is  the  product  of  the  number  of  particles  per  unit  volume,  the 
scattering  cross  section  per  particle  and  1  minus  the  average  of  the  cosine  of 
the  scattering  angle.  For  the  particles  of  interest  here  Mie  theory  gives  values 
for  this  product  of  7.43mm-l,  38.42mm-l,  42.81  mm-1,  38.50/n/n-l  and 
33.90m/w-l  for  increasing  particle  size  and  the  10%  case.  The  1%  values  are 
simply  a  10th  of  these.  From  the  figures  scatter  patterns  follow  the  general 
behaviour  of  this  parameter,  the  smallest  particle  giving  the  smallest  signal 
and  the  middle  sized  the  largest. 

Work  is  progressing  to  perform  the  experiments  for  a  wavelength  of  830nm 
(which  is  more  suitable  for  medical  applications)  and  to  obtain  calculations  of 
the  expected  scatter  pattern  using  the  diffusion  approximation. 
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Meastrements  of  angular  scattering  by  rough  surfaces  at  grazing  incidence 

V.  Ruiz  Cortdi  and  J.C.  Dainty 

Applied  Optics  Section.  Blackett  Laboratory. 

Imperial  College  of  Science,  Technology  and  Medicine 
Prince  Consort  Road,  London  SW7  2BZ,  U.K. 


An  experimental  study  of  light  scattering  at  grazing  incidence  from  random  rough  surfaces  is 
present^'4  The  surfaces  are  fabricated  in  photoresist  and  gold-coated. 


The  interaction  between  light  and  matter  has  been  studied  for  several  years.  The  reflection  and 
refraction  of  light  in  a  plane  surface  of  any  material  are  well  known  phenomena  and  it  is  possible  to  obtain 
information  about  material  proDerties  by  analysing  the  reflected  and  refracted  light.  However,  the 
interaction  of  light  with  rough  surface  is,  in  general,  a  non-resolved  problem.  In  the  last  few  years 
significant  advances  have  been  made  in  the  study  of  light  scattered  from  rough  surfaces.  Of  special  interest 
has  been  multiple  scattering  effects 1  at  normal  and  small  angles  of  incidence  (up  to  approx.  50°  from  the 
normal).  Little  work  has  been  done  at  grazing  incidence  where  the  current  theories  are  no  longer  valid. 
Furthermore,  the  available  data  at  such  angles  of  incidence  does  not  include  accurate  information  about 
the  surface  structure.  The  study  of  rough  surface  scattering  at  grazing  incidence  has  many  potential 
applications,  such  as  radar. 


Forscalarwsve  fields  at  arbitrary  angles  of  incidence  the  scattered  field  isdescribed  by  the  Helmholtz 
integral  formula2 

where  is  the  scattered  field  at  any  arbitrary  point  r„  in  the  medium,  4/(r)  is  the  incident  field  on  the 
surface,  d/dn  is  ihe  derivative  along  the  normal  component  and  G  (r0,  r )  is  the  free  space  Green ’s  function. 

Different  approximation  methods  have  been  used  to  solve  this  integral  equation  such  as  Perturbation 
Theory 3  for  small  surface  heights  and  the  Kirchhoff  Approximation 2  for  small  slopes,  each  of  these  methods 
have  a  range  of  validity  depending  on  the  parameters  of  the  rough  surface.  However  some  methods  fail 
at  large  angles  of  incidence  (greater  than  50* )  or  others,  due  to  computer  implementation,  are  excessive 
time-consuming  programs.  To  avoid  these  problems  for  grazing  incidence  it  is  possible  to  apply  the 
Parabolic  Approximation  Method 4, 
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At  grazing  incidence  the  scattered  field  propagates  predominantly  in  one  direction  (fig  1),  it  has  a 
slow  varying  component  in  the  z-axis  and  this  fact  is  used  for  the  parabolic  approximation  to  redefine  the 
Green’s  function  by  a  direct  approximation  of  the  free  space  function  which  for  a  2-D  geometry  is  given 
by4 5: 
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With  this  approximation  we  solve  the  Helmholt2  integral  equation  and  we  obtain  an  expression  for  the 
scattered  field,  of  course  this  method  is  inapplicable  to  situations  in  which  backscattering  is  significant. 
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Figure  1.-  Measured  diffuse  scattering  at  different  angles  of  incidence  (60°,  70°  and 
80°)  from  a  rough  surface  with  standard  deviation  a*  =  2.27|iwt  and  correlation 

length  a  =  20.9|im 
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A  simple  numerical  method  based  on  the  Kirchhoff  approximation  with  shadowing  is  presented  for  the  calculation 
of  the  scattered  intensity  distribution  for  rough  surfaces  illuminated  at  angles  up  to  «  83°  incidence.  It  may  be  possible 
to  go  to  higher  angles  but  man y  more  discretisation  points  are  required  on  the  surface.  The  method  is  valid  for  a  Tide 
range  of  roughness  values  which  is  an  advantage  over  other  methods. 


1.  Introduction 


2.  Theory 


Scattering  of  electromagnetic  waves  from  rough  sur¬ 
faces  for  high  angles  of  incidence  has  proved  to  be  a 
difficult  problem  to  calculate.  The  situation  is  of  in¬ 
terest  particularly  for  scattering  of  radar  from  the  sea 
surface  where  the  incidence  angle  is  within  a  few  de¬ 
grees  of  grazing  { 1  ].  Calculations  for  this  situation 
have  proved  very  difficult  although  some  have  been 
done.  The  parabolic  equation  method  [2-4]  has  been 
used  to  calculate  from  very  near  grazing  incidence  but 
it  is  limited  to  surfaces  with  height  variations  less  than 
approximately  one-fifth  of  a  wavelength  of  the  illumi¬ 
nating  radiation.  The  small  perturbation  method  [5] 
has  also  been  used  but  it  is  limited  to  even  smaller 
variations  of  the  surface  heights. 

In  this  communication  the  single  scatter  Kirchhoff 
approximation  with  shadowing  [6-1 1  ]  is  used  to  cal¬ 
culate  the  scattered  intensity  patterns  for  high  angles 
of  incidence.  The  physical  justification  for  using  this 
approach  is  that  at  high  incidence  most  of  the  sur¬ 
face  will  be  in  shadow  and  only  the  high  parts  of  the 
surface  will  be  scattering.  Since  these  will  tend  to  be 
flat  (as  they  will  be  maxima  of  the  surface)  the  light 
will  tend  to  be  scattered  away  from  the  surface.  Hence 
only  the  single  scatter  term  is  required.  The  results 
are  compared  to  experimental  results  and  the  varia¬ 
tion  of  the  scattered  intensity  with  changing  rough¬ 
ness  is  presented  to  show  that  the  method  is  valid  for 
rougher  surfaces. 


The  required  results  are  derived  from  the  two- 
dimensional  Helmholtz  integral  equation  for  the 
scattered  field  (the  surface  profile  is  constant  along 
the  y -direction) 


r 

8E, 


(1) 


where  Ei(x',z')  is  the  total  field  at  the  point  (*',  z' ) , 
r  =  [(jc-jc')2+  (z-z')2l,/J>#oU(^r)  is  the  zeroth 
order  Hankel  function  of  the  first  kind  and  ds'  is  an 
element  of  the  surface. 

The  Kirchhoff  Approximation  approximates  the  to¬ 
tal  field  on  the  surface  as  the  sum  of  the  incident  field 
plus  the  reflected  field  at  that  point 


Et(x,z)  =  (1  +  R)Ei(x,z),  (2) 

~~  =  i(l -R)k,m  Ei(x,z),  (3) 

where  R  is  the  planar  reflection  coefficient,  at  the 
point  (x,z),  which  depends  on  the  local  incidence 
angle,  E,(x,  z)  is  the  field  incident  at  that  point,  k,  is 
the  incident  wave-vector  and  it  is  the  outward  normal 
to  the  surface  at  that  point.  The  condition  for  valid¬ 
ity  of  the  approximation  is  immediately  obvious  from 
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above  since  R  is  the  planar  reflection  coefficient.  The 
surface  must  be  locally  flat  or,  as  is  usually  written, 
the  radius  of  curvature  of  the  surface  must  be  large 
compared  to  the  wavelength. 

Substituting  (2 )  and  (3)  into  ( 1 )  and  proceeding  as 
in  Beckmann  [  12]  the  standard  single  scatter  solution 
can  be  written  as 

<pt?(0i,6)  =  ^exp[i{kor-tt/4)] 

x  J cxp{i[(ko(smG,-sin0)x 
r 

-  knicosO,  +  cos#)  AU)}} 

x  [-A<)(sin0  +  sinft) h'(x) 

+  ko(cos6  -  costfi) 

-  R(.x,z)ko(smO  -sin6,)h'(x) 

+  R{x,z)ko(casd  +  cos0i)]&>c,  (4) 

where  the  reflectivity  R(x',  z‘ )  is  a  function  of  x  since 
the  local  angle  of  incidence  varies  along  the  surface. 
In  eq.  (4),  F  represents  the  surface  profile,  0,  is  the 
incidence  angle,  0  is  the  angle  of  scatter  (measured 
positive  in  the  opposite  sense  to  the  incident  angle). 
R(x,z)  is,  in  the  general  case  given  by  the  Fresnel 
reflection  coefficients.  To  obtain  this  expression  we 
have  assumed  an  incident  plane  wave  of  unit  ampli¬ 
tude.  When  the  Kirchhoff  method  is  used  this  expres¬ 
sion  is  usually  simplified  by  integration  by  parts  (6- 
9],  However  this  process  leads  to  two  terms,  one  of 
which  is  an  edge  effect  term  which  is  neglected  [6-9]. 
In  the  situations  of  very  nearly  grazing  incidence  this 
edge  term  is  not  negligible,  so  the  calculation  proceeds 
from  a  discretisation  of  the  above  equation.  This  in¬ 
volves  splitting  the  integral  into  a  summation  of  val¬ 
ues  which  are  constant  over  a  small  region  of  the  vari¬ 
able  x.  This  equation  is  the  first  term  in  the  iterative 
solution  method  [13];  however,  to  obtain  a  physically 
realistic  result  the  effects  of  shadowing  must  be  in¬ 
cluded.  This  is  particularly  true  for  the  case  of  grazing 
incidence  when  a  large  fraction  of  the  surface  is  not 
illuminated  because  the  incident  light  is  blocked  by 
other  parts  of  the  surface.  These  effects  are  included 
by  multiplying  the  integrand  in  (4)  by  incidence  and 
scatter  shadow  functions 


S(x',z')  =  1,  if  (jc',  z‘)  is  illuminated, 

=  0,  if  (x',z')  is  rot  illuminated, 

S'(x :',/)  =  1,  if  {x',z')  is  visible, 

=  0,  if  (x\  z')  is  not  visible. 

These  represent  geometrical  or  straight  line  shadow 
functions  and  so  are  an  approximation  to  the  true  ef¬ 
fect  of  shadowing.  Equation  (4)  was  discretised  by  as¬ 
suming  that  the  integrand  is  constant  over  a  suitably 
small  range  of  the  surface  profile.  Then  the  integral 
can  be  replaced  by  a  summation  and  the  system  ran 
be  programmed  into  a  computer.  To  remove  speckle 
noise  tae  calculations  are  performed  for  many  sur¬ 
faces  of  the  same  statistics  and  the  resulting  intensity 
patterns  averaged. 

It  should  be  pointed  out  that  the  inclusion  of  the 
shadowing  explicitly  is  essential  for  the  method  to 
work  for  grazing  incidence.  It  was  found  that  for  graz¬ 
ing  incidence  many  more  points  were  required  in  the 
discretisation  of  the  surface  profile  to  follow  the  phase 
variations  of  the  incident  field  on  the  surface.  This 
means  that  the  second  (and  higher)  order  terms  in 
the  iterative  series  solution  (i.e.  without  the  inclusion 
of  the  shade  wing  explicitly)  would  take  too  long  to 
compute  on  even  the  fastest  computers.  This  is  im¬ 
portant  since  it  was  shown  [13]  that,  for  the  itera¬ 
tive  solution,  if  an  incident  ray  is  blocked  by  n  other 
points  before  it  can  reach  a  particular  point  on  the 
surface  the  n’th  iteration  is  required  to  account  for 
the  shadowing.  For  example  in  fig.  1  the  ray  shown 
is  blocked  by  points  A,  B,  C  and  D  before  reaching 
point  E.  Hence  the  fourth  order  iteration  would  be  te¬ 
nured  to  give  the  correct  scattered  field  from  point  E. 
For  grazing  incidence  light  will  tend  to  be  blocked  by 
many  points  and  so  very  high  orders  in  the  iterative 
isries  would  have  to  be  computed  to  give  an  accurate 


Fig.  1.  Shadowing  at  large  angles  of  incidence.  The  incident 
ray  that  is  in  the  direction  of  point  E  is  blocked  by  the  four 
points  A,  B,  C  and  D. 
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result  Therefore,  although  the  method  presented  here 
is  an  approximation  to  the  iterative  solution,  it  is  a 
practical  solution  whereas  computer  time  limitations 
are  prohibitive  for  the  iterative  method. 


3.  Results 

In  this  section  results  are  presented  and  compared 
with  experimental  data.  The  calculations  were  per¬ 
formed  for  a  one-dimensional  perfectly  conducting 
surface  and  for  a  metal  surface  but  only  the  perfectly 
conducting  surface  results  are  presented.  For  the  metal 
case  there  was  no  difference  in  the  normalised  curves 
for  the  two  polarisations  although  the  absolute  value 
of  the  scattered  energy  was  higher  for  the  s  polarised 
case  (electric  vector  perpendicular  to  the  scattering 
plane)  than  the  p  polarised  case.  This  is  due  simply 
to  the  fact  that  the  Fresnel  reflection  coefficients  are 
larger  for  s.  The  normalised  curves  were  very  similar 
to  the  perfect  conductor  case.  For  a  perfect  conductor 
the  scattered  intensity  is  the  same  for  the  two  funda¬ 
mental  polarisations,  j  and  p,  and  the  modulus  of  the 
reflection  coefficient  is  1. 

The  surface  used  here  is  approximately  described 
by  a  gauss  ian  correlation  function  with  a  1  /e  distance 
t  as  12.0  pm  and  a  gaussian  distribution  of  heights 
with  a  standard  deviation  of  <r  =  0.40  pm.  The  ex¬ 
periments  were  conducted  with  a  helium-neon  laser 
(X  =  0.633  /rm).  The  calculations  were  performed  on 
a  4002  length  of  surface  discretised  into  1000  points. 
On  a  Sun  Sparc  station  the  calculations  took  40  min¬ 
utes  per  frame.  As  a  check  on  the  validity  of  the  cal¬ 
culation  the  ratio  of  the  scattered  energy  and  the  in¬ 
cident  energy  was  found  This  is  termed  the  unitarity 
as  it  should  have  the  value  1  for  a  perfect  calculation 
(as  all  the  energy  is  reflected  from  a  perfectly  conduct¬ 
ing  surface).  The  unitarities  were:  for  60°  incidenc- 
=  1.0089;  for  70°  =  0.9978  and  for  80°  =  1.0436. 
These  numbers  compare  well  with  results  for  other 
situations  [6-8}  and  give  confidence  in  the  shape  of 
the  resulting  curves. 

Figure  2  shows  the  comparison  between  calcula¬ 
tions  and  experimental  results  for  the  normalised  in¬ 
tensity  (the  graphs  are  normalized  to  have  unit  area 
under  the  curve).  Note  that  the  experimental  results 
are  for  a  surface  rough  in  both  the  x  and  y  directions 
whereas  the  calculation  is  for  a  surface  rough  in  the  x 


scatter  angle  (degrees) 


Fig.  2.  Comparison  of  experiments  and  calculations  for  60°, 
70°  and  80°  incidence  on  a  surface  with  <j  «=  0.4  pm  and 
t  “  12.0  pm.  Solid  lines  are  the  calculation  and  circles 
are  the  experimentally  measured  values  for  p  polarisation 
incident  and  detected.  All  curves  are  normalised  to  unit  area. 

direction  and  constant  i.’  the  y  direction.  This  means 
that  the  quantitative  comparison  of  the  results  may 
not  be  valid.  However,  the  qualitative  comparison  can 
be  useful.  As  can  be  seen  from  the  figure  these  curves 
for  the  scattered  light  distributions  as  functions  of  an¬ 
gle  agree  well  with  the  experimental  data.  In  particu¬ 
lar  the  presence  of  the  large  specular  peak  at  80°  in¬ 
cidence  is  predicted  in  the  calculations. 

For  higher  angles  of  incidence  the  calculations  the 
method  starts  to  break  down.  The  values  for  the  uni¬ 
tarity  are  1.02  for  85°  incidence,  1.06  for  86°,  1.20  for 
87°  and  2.51  for  89°.  Therefore  for  the  values  of  the 
surface  parameters  used  the  met'  oti  described  here 
may  be  used  for  incidence  angles  of  up  to  approxi¬ 
mately  86°.  It  may  be  that  very  many  extra  points  are 
required  for  the  large  angles  of  incidence  to  follow  the 
phase  of  the  incident  wave  on  the  surface. 

In  the  computer  it  is  possible  to  increase  the  rough¬ 
ness  of  the  surface  by  increasing  the  value  of  o.  Figure 
3  shows  the  normalised  scattered  intensity  patterns  for 
four  surfaces  all  with  the  same  value  of  r  =  12.0  lira 
as  for  the  surface  above  and  values  of  the  roughness  of 
(a)  u  =  0.40  pm,  (b)  c  =  0  8  pm,  (c)  a  =  1.2  pm 
and  (d)  a  =  1.6  pm.  The  incidence  angle  is  80°  for 
the  curves  shown.  The  unitarities  for  (he  calculations 
were  (a)  1.0436,  (b)  0.96,  (c)  0.998  and  (d)  1.02. 
All  of  these  numbers  are  reasonable  and  show  that  the 
calculation  is  less  sensitive  to  the  roughness  than  it  is 
to  the  angle  of  incidence  when  the  incidence  angle  is 
high.  This  is  true  since  most  of  the  surface  is  shadowed 
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scatter  angle  (degree!) 


Fig.  3.  Calculated  variation  of  scattering  for  increasing 
roughness.  The  angle  of  incidence  is  80°.  The  solid  line  is 
for  a  *  0.4  ftm,  the  dashed  line  for  a  -  0.8  /im,  the  crosses 
for  a  =  1.2  ftm  and  the  circles  for  o  =  1.6  /an,  all  for 
r  =  12.0  /an. 


is  valid  for  incidence  angles  up  to  «  85°  and  for  sur¬ 
faces  whose  rms  roughness  is  up  to  a  few  wavelengths 
and  can  be  used  to  calculate  the  scattered  intensity 
for  any  material  of  rough  surface.  The  method  has  the 
virtue  of  simplicity  and  seems  to  agree  with  experi¬ 
mental  data. 
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